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The conditions for synchronization in unidirectionally coupled chaotic oscillators are revisited. 
We demonstrate with typical examples that the conditional Lyapunov exponents (CLEs) play an 
important role in distinguishing between intermittent and permanent synchronizations, when the 
analytic conditions for chaos synchronization are not uniformly obeyed. We show that intermittent 
synchronization can occur when CLEs are very small positive or negative values close to zero while 
permanent synchronization occurs when CLEs take sufficiently large negative values. There is also 
strong evidence for the fact that for permanent synchronization the time of synchronization is 
relatively low while it is high for intermittent synchronization. 
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I. INTRODUCTION 



Much attention has been focussed on the synchro- 
nization of chaotic systems through different coupling 
schemes during the past decade or so. It has potential 
applications not only in secure communication but also 
in chaotic cryptography]!]]. In order to achieve synchro- 
nization, coupled chaotic systems have to satisfy certain 
conditions. According to Pecora and Carroll the cou- 
pled chaotic systems can be regarded as master and slave 
systems which will perfectly synchronize only if the sub- 
Lyapunov or conditional Lyapunov exponents (CLEs) are 
all negative J|, [|. However, recently there are reports 
claiming that synchronization can be achieved even with 
positive CLEs||, |[ Also there are situations where 
one can observe intermittent or round-off induced syn- 
chronization phenomenon |t], |§|, P where locking occurs 
at certain time intervals only, so that the synchroniza- 
tion of chaos need not be always permanent. It is there- 
fore important to understand clearly the conditions un- 
der which chaos synchronization can occur and to know 
how to distinguish between permanent and intermittent 
synchronizations exhibited by coupled chaotic systems. 

In this paper, the analytic conditions for synchroniza- 
tion in coupled chaotic systems are revisited to show 
the difficulties in distinguishing between intermittent and 
permanent synchronizations which we illustrate by using 
simple chaotic systems. We demonstrate that the con- 
ditional Lyapunov exponents (CLEs) play an important 
role in distinguishing between intermittent and perma- 
nent synchronization, when the analytic conditions for 
chaos synchronization are not uniformly obeyed. We also 
point out that intermittent synchronization can occur 
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when the largest CLE has a value close to zero (either 
positive or negative) while permanent synchronization 
occurs for sufficiently large negative values of the CLEs. 

In Sec. II, we describe briefly the notion of chaos syn- 
chronization and the analytic conditions for synchroniza- 
tion which are to be uniformly obeyed as well as the 
existence of negative conditional Lyapunov exponents, 
with reference to a simple coupled system of autonomous 
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Duffing- van der Pol (ADVP) oscillators JToJ, |TT). Sec 
describes how the CLEs can be used to analyse the na- 
ture of synchronization in coupled ADVP oscillators. In 
Sec. we show how conditional Lyapunov exponents 
play an important role in distinguishing between the 
intermittent and permanent synchronizations occurring 
in a class of unidirectionally coupled nonlinear oscilla- 
tors including coupled chaotic pendula|@], coupled Duff- 
ing oscillators ra and coupled Murali-Lakshmanan-Chua 
(MLC) circuitspO[ either in the absence of the analytic 
criteria or when they are not uniformly obeyed. We also 
show how the time of synchronization can be used quali- 
tatively to check the presence or absence of intermittent 
synchronization. Sec. ^ summarizes the results. 



II. CHAOS SYNCHRONIZATION AND THE 
CONDITIONS 



In this section, we consider a set of two unidirection- 
ally coupled autonomous Duffing-van der Pol (ADVP) 
oscillators as a model system to analyse the concept of 
chaos synchronization. Further, we show that the derived 
analytic conditions for the occurrence of chaos synchro- 
nization are uniformly obeyed in this system for a specific 
set of parametric values (that is, the analytic conditions 
are shown to be valid for all values of the variables, say, 
(xi,X2,xs) € 1Z 3 for all t). We also indicate that when 
the usual analytic conditions are not uniformly obeyed by 
a given system, then it may either exhibit intermittent 
or permanent synchronization^. In this case, the na- 
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ture of the conditional Lyapunov exponents can be used 
profitably to distinguish between intermittent and per- 
manent synchronization. 




100 



FIG. 1: Figure showing the variation of synchronization error 
rj as time elapses in ADVP oscillators [Eqs. ([!])] for v — 100, 
a = 0.35, = 300 and e = 1.35. 

We consider the ADVP model in which Murali and 
Lakshmanan demonstrated chaos synchronization ]Tc| , 
|TT[ . The rescaled and dimensionless version of the unidi- 
rectionally coupled ADVP oscillators can be written as 
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Here v, a, and /3 are rescaled parameters, which are fixed 
at v = 100, a — 0.35 and j3 — 300. One can define a 
measure of the synchronization error, 77, as 



V = \J {x m - Xs) 2 + {y m ~ Vs) 2 + (z m - Z s f 



(2) 



For synchronization, the above measure r\ —> as t — » 
00. The synchronization error rj versus time is plotted 
in Fig. [l| The falling up of r\ to zero is an indication of 
synchronization at a finite time. However, this alone does 
not ensure that the synchronization is permanent and 
one has to verify additional criteria. For this purpose, 
we will first discuss conditions for which one can achieve 
permanent synchronization. 

(i). The criterion introduced by Fujisaka and Yamada 
fl2[ 13 1 for high quality synchronization requires 
that the largest eigenvalue of the Jacobian matrix 
corresponding to the flow evaluated on the synchro- 
nization manifold be negative. In order to check 
this for the system (|l]), let us consider the spe- 
cific choice e = 1 + a = 1.35 in Eqs. (pi). In 
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this case, then one can write the difference sys- 
tem of the ADVP oscillators for x* = (x m — x s ), 
y* = (y m — Vs), z* = (z m — z a ) in matrix form as 



(3) 



where a — (x m — x s ) 2 + 3x m x s > 0. The slave 
system l^b) synchronizes perfectly with the master 
system ([Hi ) only if all the eigenvalues of the above 
linear system (^) possess negative real parts. One 
can easily prove that this is indeed the case for (||). 

Recently, He and Vaidya developed a crite- 
rion for chaos synchronization based on the no- 
tion of asymptotic stability of dynamical systems, 
which refers to the condition for a given chaotic 
system with master-slave configuration to reach the 
same eventual state at a fixed (but sufficiently far 
enough) time irrespective of the choice of initial 
conditions. One of the practical ways to establish 
the asymptotic stability of the response subsystem 
is to find a suitable Lyapunov function which can be 
defined as the square of the magnitude of the vector 
describing the distance from the synchronization 
manifold. Then the condition for all the pertur- 
bations to decay to the synchronization manifold, 
without transient growth, is that the time rate of 
the Lyapunov function has a negative magnitude 
for all times p5|. 

Now for the difference system (^), one can analyt- 
ically write a Lyapunov function in the following 
form for (3,v> 0@ [tl|, 



E 



vy 



*21 



,*2 



> 0. 



(4) 



Then the rate of change of E along the trajectories 
is given by 



E 



uf3 [ax* 2 + (x* -y*) 2 ] < 0. 



(5) 



Since E is a positive definite function and E is neg- 
ative definite for (x*,y*,z*) e 1Z 3 , according to 
Lyapunov theorem x* , y* and z* — > as t — > 00. 
Thus perfect synchronization occurs as t — > 00. 

In addition, if all the sub-Lyapunov or conditional 
Lyapunov exponents (CLE) are negative then one 
may have perfect synchronization between the mas- 
ter and slave systems. The CLEs are the corre- 
sponding Lyapunov exponents of the slave system. 
For the coupled ADVP oscillators [Eqs. (Q)] the 
CLEs can be calculated numerically using the stan- 
dard algorithm [fl6|. We find that for the system (|l|) 
with e = 1 + a = 1.35, the numerical value of the 
largest CLE is —0.3134. Thus, one may expect that 
the ADVP oscillators will synchronize perfectly for 
e = 1.35 which is indeed true. 
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Consequently the conditions for synchronization in cou- 
pled chaotic systems may be atleast any one of the fol- 
lowing: 

(i) . The largest eigenvalue of the Jacobian matrix corre- 

sponding to the flow evaluated on the synchroniza- 
tion manifold must always be negative [see Eq. @)]. 

(ii) . Existence of a suitable Lyapunov function for the 

difference system as discussed earlier. 

(iii) . The sub- Lyapunov exponents or CLEs are all neg- 

ative. 

Among these three condi tion s, it is obvious that either 
of the conditions |ij and (ii] is both necessary and suf- 
ficient for perfect or permanent synchronization as the 
very definition of the latter implies these conditions. On 
the other hand the condition (iii) is only a necessary one, 
because the CLEs pertain to finite time averages only. 
As we have discussed above for the specific parametric 
(e = 1.35) case of the coupled AD VP oscillators (|l|), all 
the conditions (i) - (iii) are satisfied and so pemanent 
synchronization indeed occurs. However, suppose for a 
given system, if the first two conditions are not uniformly 
obeyed and only the third condition is satisfied, the ques- 
tion is whether the negativity of CLEs alone ensures per- 
manent synchronization. 

In fact, recent studies show that in certain specific sys- 
tems the observed synchronization is not always a per- 
fect one and rather one can have an intermittent syn- 
chronization^ [TtJ even when all the CLEs are negative. 
Further, there are reports stating that one can achieve 
synchronization in the presence of positive conditional 
Lyapunov exponents as well[Q. So one has to clarify 
the conditions under which permanent synchronization 
can occur in coupled chaotic systems and how one can 
distinguish between the intermittent and permanent syn- 
chronization. 



III. NATURE OF CLES IN ADVP 
OSCILLATORS [Eqs. (Q)] FOR GENERAL e 

Based on the analysis of the previous section, one un- 
derstands that the coupled ADVP oscillators (0) show 
perfect synchronization for e = 1.35. Now what hap- 
pens to other values of e in the coupled system (|])? It 
does not seem to be feasible either to deduce analyti- 
cally the nature of the largest eigenvalue of the Jacobian 
corresponding to the flow evaluated on the synchroniza- 
tion manifold or to obtain explicitly a suitable Lyapunov 
function. 

In order to understand the nature of synchronization 
in Eqs. (|l|) for e ^ 1.35, let us calcuate the CLEs. Fig. ^ 
shows the variation of the largest CLE as a function of the 
coupling strength (e). From the figure, one may expect 
that the ADVP oscillators will synchronize for e > 0.65. 
However, the largest CLE value is almost zero in the 




FIG. 2: Figure showing the variation of largest CLE as a 
function of coupling strength (s) for ADVP oscillators. 



range of e <E (0.5,0.65). On careful observation, by ex- 
amining the synchronization error r\ with the addition of 
a small amount of noise at every integration step, one 
finds that perfect synchronization will occur in the cou- 
pled ADVP oscillators for e > 0.65 only. 

In addition, the Jacobian matrix corresponding to the 
flow evaluated on the synchronization manifold can be 
given by 



-"(34n + c) v 

J= I 1 -1 -1 

p 



(6) 



where e = a + c (c = constant), 
eigenvalues take the form 



a = 0.35. Then the 



Al,: 



1 + v(3x 2 sm + c) 



± 



vTi + K3z s 2 m + c)] 2 - M3ZL + c - 1) 



A, = 0. 



(7) 



In the above, the sign of the largest eigenvaule is de- 
termined by the factor (3Xg m + c — 1) inside the square 
root sign. When c = 1 (that is, e — 1.35), the largest 
eigenvalue is negative as discussed earlier. Further, it is 
evident that the largest eigenvalue remains negative for 
c> 1. 

On the other hand, if c < 1 then one has to eval- 
uate the eigenvalues numerically in order to check the 
nature of synchronization. For this purpose we exam- 
ine the variation of the eigenvalues numerically. Figs. ^ 
show the variation of the largest eigenvalue as a function 
of time. Fig. |](a) illustrates the variation of the largest 
eigenvalue for e — 0.60 (c — 0.25). Here the largest eigen- 
value oscillates chaotically about zero with an average of 
the order of +10 -1 and hence no synchronization is pos- 
sible for e — 0.60. For e — 0.65, the largest eigenvalue 
still oscillates chaotically about zero [see Fig. |(b)] with 
an average « 10~ 2 . In the latter case the largest eigen- 
value spends more time in the negative region despite 
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function corresponding to the square of the magnitude of 
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FIG. 3: Figure showing the plot of the largest eigenvalue of 
the Jacobian matrix (^|) versus time for AD VP oscillators: (a) 
e = 0.60, (b) e = 0.65, (c) e = 1.00 and (d) e = 1.35. 



the fact that its average value is positive. In addition, 
for e = 0.65, the CLE has a value of -0.01298 and thus 
it shows permanent synchronization. Similar arguments 
hold good for 0.65 < e < 1.35. Figs. |(c) and |(d) de- 
pict the largest eigenvalue for e = 1.00 and e = 1.35, 
respectively. 

In addition to the above, one can also verify the con- 
dition (ii) , by evaluating the time rate of the Lyapunov 
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FIG. 4: Figure showing the plot of E in the AD VP oscillators 
[Eq. | asa function of time for (a) e = 0.64 and (b) e = 0.65. 



the vector describing the distance from the synchroniza- 
tion manifold |15|, which is given by 



*2 



E = -v(3x sm + c)x*~ - y 

+{l + v)x*y* + (p + l)y*z* <0. 



(8) 



Fig. |4|(a) shows the variation of E for e = 0.64. From 



the figure it is clear that the condition [ii] is not uni- 
formly satisfied. Thus, there is no synchronization in the 
coupled AD VP oscillators ([[]) for e — 0.64. However, 
E — » at a finite time for e = 0.65 and e > 0.65 [see 
Fig. f|(b)] and hence perfect synchronization does occur 
in the coupled AD VP oscillators for e > 0.65. 

Thus from the above analysis and from the nature of 
the CLEs, one may conclude that perfect synchronization 
in the coupled ADVP oscillators (|l|) will occur for e > 
0.65. The role of CLEs become even mor e im portant 
when both the analytic conditions |(ij and (ii) arc not 
uniformly satisfied. In such cases, in the following, we 
will analyse how the nature of CLEs can be effectively 
used to understand synchronization in typical coupled 
chaotic systems. 
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IV. DISTINGUISHING INTERMITTENT AND 
PERMANENT SYNCHRONIZATION USING 
CLES 

We now wish to clarify the conditions for which one can 
observe intermittent and permanent synchronization, al- 
beit in a qualitative way, when the analytic conditions |^| 
and (ii) are not uniformly obeyed. In this regard we wish 
to analyse the nature of synchronization on the basis of 
the conditional Lyapunov exponents (CLEs) by consid- 
ering three typical dynamical systems as examples. They 
are namely (i) coupled chaotic pendula, (ii) coupled Duff- 
ing oscillators and (iii) coupled MLC circuits. We find 
that all these thre e sy stems do not satisfy the analytic 
conditions (i) and (ii) uniformly in contrast to the cou- 
pled AD VP oscillators considered in Sec. II. We demon- 
strate how the nature of the CLEs become important in 
distinguishing the intermittent and permanent synchro- 
nizations exhibited by these systems. 



A. Coupled chaotic pendula 

First we consider a pair of coupled chaotic pendula 
defined by the following set of equations and investigated 
by Baker, Blackburn and Smithj^] 



Tq cos Qi, 

- 1 e a = r cosftt + 



c sin a . — sin ( 



(9a) 
(9b) 



Here 9 m and 9 S correspond to the angular coordinates of 
the master and slave systems, respectively. Q~ l corre- 
sponds to the damping factor, Tq is the normalized drive 
torque, Q is the drive frequency and c is the coupling 
strength. By fixing the parameters at Q = 5.0, To = 1.2 
and Q = 0.5, the quality of synchronization in the above 
coupled system @ has been studied. This has been fa- 
cilitated by finding the synchronization error, 



1 1 



(10) 



oo. 



As seen earlier, for synchronization rj — > as t 

It has been shown by Baker, Blackburn and Smith|7j 
that permanent synchronization in the above coupled 
identical chaotic pendula (||) does not occur except as a 
numerical artifact arising from finite computational pre- 
cision. Further, they showed that the synchronization 
of the above coupled pendula ([|) is always intermittent, 
for any value of the coupling coefficient, by using nu- 
merical and analytical tests. However, contrary to the 
above the present authors have reported that there ex- 
ists atleast certain range of coupling coefficient values 
for which one can observe permanent synchronization by 
computing the CLEs[^7[ [ll|. In order to understand this 
we proceed as follows. 

For the coupled chaotic pendula ([)]), as shown in 
Ref. R| the eigenvalues of the Jacobian matrix corre- 
sponding to the flow evaluated on the synchronization 



manifold take the form 
-1 

2 ~ 2Q 
A 3 = 0. 



Ai. 



1 ± y/l-2Q 2 (l-c) cos9 s 



(11) 



Here 9 sm corresponds to the angle coordinate on the syn- 
chronization manifold. On careful numerical observation, 
it has been shown by Baker, Blackburn and Smith j?J 
that the term inside the square root sign varies chaoti- 
cally about unity with an average value which is slightly 
less than one for a range of c values. This implies that 
the largest eigenvalue is not always less than zero and 
hence the condition |ij is not obeyed uniformly. 

Further, calculating the Lyapunov function which is 
the square of the magnitude of the vector describing the 
distance from the synchronization manifold fll5|, the suf- 
ficient condition for permanent synchronization can be 
written as[Q 

E = (1 - c - cos sm ) 9*uo* - Q^uj* 2 < 0, (12) 

where u = 0, 9* = 9 m — 9 S and lu* — u> m — u) a . The time 
series of the above expression varies intermittently about 
zero with an average slightly less than zero for a ran ge 
of c values [c G (0.79, 1.0)] and hence the condition [ii] is 
also not obeyed for the coupled chaotic pendula. Thus, 
in the light of the above evidence, Baker, Blackburn and 
Smith iQ, H conclude that intermittent synchronization 
can be a plausible behaviour in the coupled pendulum as 
the locking occurs intermittently. In order to verify this 
assertion, we have carried out an analysis jl7] by calculat- 
ing the conditional Lyapunov exponents of system (^J) in 
the parameter range c 6 (0.75,1.0). We find that there 
exists a range of c values for which one can indeed have 
permanent synchronization, where Baker, Blackburn and 
Smith have expected intermittent synchronization. In 
our calculations we have used the standard Wolf et al 
algorithm JT(| and used 5000 drive cycles for calculations 
after leaving out 5000 drive cycles as transient. 

Figure ||(a) shows the variation of CLE as a function 
of c for the coupled pendula and it takes negative values 
for 0.796 < c < 1.0 only. The CLE value for c = 0.79 
is +0.00320 (« 10~ 3 ), for which synchronization can not 
occur and hence the observed intermittent synchroniza- 
tion (see Fig. 1 in Ref |Q) is a computer artifact. How- 
ever, we find that for c = 0.932, the CLE value becomes 
the lowest (—0.06825) and it is relatively a large negative 
value for which intermittent synchronization is absent as 
shown in Fig. [f|(b) . The absence of the intermittent syn- 
chronization has been verified upto 2 x 10 6 drive cycles 
even with the addition of tiny noise levels showing that 
permanent synchronization does occur for this value of c. 
The same phenomenon persists over a range of c values 
close to 0.932. Figs. |^ show the variation of the synchro- 
nization error for various c values. From the figures it 
is clear that one can indeed have intermittent synchro- 
nization for a range of c values as noted by Baker, Black- 
burn and Smith[[7|. For example, the system definitely 
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FIG. 5: (a) Variation of the conditional Lyapunov exponent (CLE) as a function of c and (b) synchronization error rj [see 
Eq. (|10[)] versus time in drive cycles for coupled pendula. 
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FIG. 6: Figure showing the variation of synchronization error rj versus time in the coupled chaotic pendula for (a) c = 0.80, 
(b) c = 0.81, (c) c = 0.89, (d) c = 0.92, (e) c = 0.94 and (f) c = 0.98. Fluctuation in rj shows the possibility of having 
intermittent synchronization. In calculating rj, tiny noise levels (« 10 -15 ) are included at each time step to avoid round off 
induced synchronization. 



exhibits intermittent synchronization for c = 0.80, 0.81 
and 0.98, where we find that in this range of c values (that 
is, c < 0.88 or c > 0.95) the largest CLEs take values very 
close to zero (either positive or negative). But more in- 
terestingly we also find that there exists another range 
of c values from 0.89 to 0.94 where one can haveperma- 
nent synchronization. In calculating rj for Figs. o(a)-(f), 
we have added a small amount of noise (w 10~ 15 ) at 
each time step to avoid the effect of round off induced 
synchronization. We note that the CLEs are negative 
and the largest CLE is relatively away from zero in this 
range. Thus we note that the actual values of the CLEs, 
including the largest CLE, seem to distinguish between 
permanent and intermittent synchronization. 

Next we calculate the time of synchronization (t s ), that 
is, the time taken by the system for which the synchro- 
nization error becomes zero. Fig. [t] depicts the variation 
of t s as a function of the coupling coefficient. It can be 



easily seen that the value of t s is very low for the range 
of c values where the CLEs are large negative which cor- 
responds to permanent synchronization. Similarly, the 
value of t s becomes high in the regions where the largest 
CLE has a value close to zero (either negative or posi- 
tive) which corresponds to intermittent synchronization, 
as there exists large fluctuations in the finite time Lya- 
punov exponents here. 

It is clear from the above analysis on the coupled 
chaotic pendula (^) that there exists intermittent syn- 
chronization for certain ranges of coupling strength val- 
ues. However, permanent synchronization does occur for 
certain other range of c values. In particular, we have 
pointed out that when the CLEs become very close to 
zero (cither positive or negative) one can have intermit- 
tent synchronization while permanent synchronization 
does occur when the CLEs become large negative. Thus 
the CLEs play a crucial role in distinguishing between 
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FIG. 7: Figure showing the variation of time of synchroniza- 
tion t s (in drive cycles) versus coupling strength for coupled 
chaotic pendula. 



intermittent and permanent synchronization. In the fol- 
lowing subsections we will examine the existence of sim- 
ilar intermittent synchronization in the coupled Duffing 
equations and coupled MLC circuits. 



B. Coupled Duffing oscillators 

Now we consider the case of two coupled Duffing oscil- 
lators described by the following set of equations, 

x m + ax m + S3x 3 m = / cos wt, (13a) 
x s + ax s + [3x 5 s = f cosojt + c(x s — x m ), (13b) 

where a — 0.1, (3 = 1.0, / = 10 and u = 1.0 and c is 
the coupling strength. By setting x = y, z — cut, one 
can rewrite the above equations as a set of first order 
ordinary differential equations of the form 



x Ti 

Vr, 



Vs 



= Vm 

= -ayn 



fix* + / COS Z r , 



-ay s - (3x1 + f cos z s + c(x s 



(14a) 



(14b) 



First, we study the nature of synchronization in the 
above coupled Duffi ng equations by analysing the con- 
ditions (i) and (ii) of Sec. II as in the case of chaotic 
pendula discussed above. Here the eigenvalues of the Ja- 
cobian matrix corresponding to the flow evaluated on the 
synchronization manifold can be explicitly written as 



Al,2 — 



: ± y/a* + 4 (3/?a:g m + c) 
2a 
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FIG. 8: (a) Variation of the conditional Lyapunov exponent 
(CLE) as a function of c and (b) synchronization error r\ for 
coupled Duffing oscillators. 



where x sm denotes the ^-variable evaluated on the syn- 
chonization manifold. We found that the largest of eigen- 
value oscillates chaotically about zero with an average 
value very close to zero (w ±10~ 3 ) for two ranges of 
c e (1.48, 1.64) and c G (2.74,3.84). Thus the condition 
[(!)[ is not uniformly obeyed. In addition, one can write 
the condition (ii) as 



E = (l- 3/3< m - c)x*y* - ay* 2 < 0, 



(16) 



A, = 0, 



(15) 



with x* = x m —x s and y* — y m —y s . On actual numerical 
simulation, it has been found that the time series of the 
above expression oscillates chaotically about zero with 
an average of the order of ±10~ 3 . This confirms that 
the condition |(ii)| is also not satisfied uniformly. Thus, 
the only way to understand the nature of synchroniza- 
tion in Eqs. (jlj) is to analyse the CLEs. We have again 
calculated the largest CLE as a function of the coupling 
strength. Fig. || shows a plot of the largest CLE versus 
coupling strength (c). We find that the CLE takes small 
negative values close to zero (« 10 -3 ) for c € (1-48, 1.64), 
where synchronization is essentially intermittent and nu- 
merical artifact does arise in this range. In addition, 
there exists another range of c values, 2.74 < c < 3.84, 
where the CLEs are again negative (cf. Fig. ||(a)). We 
find that for c = 3.06, the CLE takes the lowest value of 
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FIG. 9: Figure showing the variation of synchronization error r\ versus time in the coupled Duffing oscillators for (a) c = 2.74, 
(b) c = 2.92, (c) c = 3.36 and (d) c = 3.84 with the inclusion of tiny noise levels of the order of 10~ 15 at each integration step. 
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FIG. 10: Figures showing the variation of time of synchro- 
nization i s versus c in the coupled Duffing oscillators. 



—0.05967 (large negative value) for which persistent syn- 
chronization occurs and the intermittent synchronization 
is absent (cf. Fig. 0(b)). The variation of synchroniza- 
tion error for a range of selected c values is depicted in 
Figs. I§](a)-(d). It can be noted from the figures that for 
c = 2.74 and 3.84 (Fig. |(a) and (d)) the syste m can 
exhibit intermittent synchronization while for c = 2.92 
and 3.36 (Fig. ^(b) and (c)) the system exhibits perma- 
nent synchronization. Here also tiny noise levels in all 
the variables were included at each integration step to 



ensure the absence of round off induced synchronization. 

As in the case of coupled pendula, we have calculated 
the time of synchronization (t s ) by varying the coupling 
strength. Fig. [h] shows the variation of t s as a function 
of the coupling strength for this case. From the figure, it 
can be noted again that t s takes relatively low values for 
a region of c values where the CLEs have large negative 
values which corresponds to permanent synchronization. 
Similarly, the values of t s becomes high in the regions 
where the largest CLE has a value close to zero or positive 
which corresponds to intermittent synchronization. 



C. Coupled MLC circuits 

Finally we consider the case of the coupled Murali- 
Lakshmanan-Chua (MLC) circuits, represented by the 
following set of equations fllOfl 



x n 



Vs 



Urn ^(^m)} 

-f3x m + cry,, 
to 



F sin z r , 



y s - h{x s ) + e(x m - x s ), 
—(3x s + ay s + F sin z s 



(17a) 



(17b) 



where e is the coupling strength and h{x) = bx — 0.5(a — 
6)(|x + 1| — \x — 1|). The uncoupled system of the above 
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FIG. 11: Figure showing (a) the variation of largest CLE as a function of the coupling coefficient e for the coupled MLC circuits 
and (b) blow up of a small region in (a) indicating the fluctuations of the largest CLE near zero. 
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FIG. 12: Figure showing the variation of synchronization error r\ versus time in the coupled MLC circuits for (a) e = 0.0377, 
(b) e = 0.0380, (c) e = 0.0390, (d) e = 0.0400, (e) e = 0.0500 and (f) a small region in (e) with the inclusion of tiny noise levels 
at each integration step. 



MLC circuits has been shown to exhibit chaos for the 
choice of parameter values j3 = 1, a = 1.015, a = —1.02 
and b = —0.55. We are interested to analyse the existence 
of intermittent synchronization in this case also. 



First let us check the validity of the conditions (i) and 



In this case, the eigenvalues of the Jacobian matrix 
of the flow evaluated on the synchronization manifold can 
be written as 

Ai, 2 = ~ [h'(x sm ) + p - e) 



±-\f[h'(x sm )+P-eY 



A, = 0. 



■A[(h'-e)0 + a], 

(18) 



Here h'(x sm ) corresponds to the derivative of h(x) with 
respect to x variable evaluated on the synchronization 
manifold. From numerical analysis, we have found that 
the largest of the above eigenvalues oscillates chaotically 
about zero with an average slightly less than zero for a 
range of e 6 (0.036, 0.04) and thus the condition Fil is 



not uniformly obeyed. On the other hand, the condition 
(ii) for the coupled MLC circuits can be written as 



E = [s - ti(x sm )} x* 2 ~ f3y* 2 + (1 - *)x*y* < 0, (19) 



with x* — x m — x s and y* = y m — y s . Again by nu- 
merical analysis, we find that the above inequality is not 
uniformly obeyed and hence the condition (ii) is not sat- 
isfied. However, the coupled syste m (|l7| ) has been shown 
to exhibit perfect synchronization[10|. In order to anal- 
yse this we have evaluated the CLEs as in the case of 
chaotic pendula and Duffing oscillators. 

The variation of the CLEs as a function of the cou- 
pling strength (e) is shown in Fig. 111]. The largest 
CLE changes its sign from positive to negative values at 
e = 0.0377 and it is slightly less than zero over the range 
0.0377 < e < 0.04. In this range the synchronization 
completely depends on the choice of initial conditions 
and the computer precision used. This implies that the 
observed synchronization is essentially due to the finite 
computational precision and accordingly locking (77 = 0) 
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occurs after relatively long time, leading to an intermit- 
tent synchronization in this range of coupling strength. 

We have also calculated the synchronization error (77) 
for various values of coupling strength as a function of 
time as in the case of previous systems and the results 
are depicted in Fig. [ljf From Figs. |l2](a)-(d), it is evi- 
dent that one can have intermittent synchronization for 
e = 0.0377, 0.0380, 0.0390 and 0.04. Fig. ||(e) and 
Fig. |l2|(f) indicates the presence of permanent synchro- 
nization for e — 0.05 where the sharp and quick fall 
off in the synchronization error towards zero occurs (see 
Fig. |H(f)). 
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FIG. 13: Figure showing the variation of time of synchroniza- 
tion t s versus coupling strength e for coupled MLC circuits. 

Fig. [l^ shows the variation of time of synchronization 
(t s ) as a function of the coupling strength (e). As in the 
case of coupled pendula and Duffing oscillators, t s takes 
relatively low values for the range of e values where the 
CLEs take large negative values corresponding to per- 
manent synchronization. However, the values of t s be- 
comes high in the regions where the largest CLE has a 
value close to zero or positive which corresponds to in- 
termittent synchronization. Thus one can conclude that 




permanent synchronization in coupled MLC circuits (|1 
occur for e > 0.04. 



V. SUMMARY AND CONCLUSIONS 

By considering typical examples of coupled oscillator 
models, we have discussed the quality of synchronization 
using various analytical and numerical tests. It has been 
noted that one can have an intermittent as well as per- 
manent synchronization in the coupled oscillator systems 
depending on the choice of coupling strength. We have 
pointed out that the Conditional Lyapunov Exponents 
(CLEs) can be used to distinguish between the intermit- 
tent and permanent synchronization when the other cri- 
teria for asymptotic stability are not uniformly obeyed. 
Particularly, we find that intermittent synchronization 
can occur when CLEs are very small positive or negative 
values close to zero while persistent synchronization oc- 
curs when CLEs take sufficiently large negative values. 
This fact is further supported by the relative time taken 
for the system to approach the synchronization manifold. 

The present work is mainly concerned with the qualita- 
tive analysis of chaos synchronization in unidirectionally 
coupled systems based on conditional Lyapunov expo- 
nents. In order to understand the entire dynamics one 
has to analyse the nature of the attractors that exist 
in the coupled system and their bifurcations. It is ob- 
vious that we have not made any quantitative analysis 
on how small or large be the magnitude of the largest 
conditional Lyapunov exponent in distinguishing the in- 
termittent and permanent synchronizations. Work is in 
progress along these directions. 
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